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ABSTRACT 



In a hot accretion flow, the radiation from the innermost region of the flow 
propagates outward and heats the electrons at large radii via Compton scattering. 
It has been shown in previous works that if the radiation is strong enough, L > 
Edd, the electrons at the Bondi radius (r B ~ 10 5 r s ) will be heated to be 
above the virial temperature thus the accretion will be stopped. The accretion 
will recover after the gas cools down. This results in the oscillation of the black 
hole activity. In this paper we show that this mechanism is the origin of the 
intermittent activity of some compact young radio sources. Such intermittency is 
required to explain the population of these sources. We calculate the timescales of 
the black hole oscillation and find that the durations of active and inactive phases 
are 3 x lO 4 (O.l/a)(M/lO 8 M )(L/2%L Edd )" 1/2 yr and 10 5 (a/0.1)(M/10 8 M o ) yr, 
respectively, consistent with those required to explain observations. Such kind 
^ ■ of feedback occurring at parsec scale should be common in low-luminosity AGNs 

and should be considered when we consider their matter and energy output. 

Subject headings: accretion, accretion discs - black hole physics - galaxies: active 

• rH ' 
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1. Introduction: intermittent activity of compact radio sources 

Feedback from active galactic nuclei (AGNs) is now widely believed to play an important 
role in the formation and evolution of galaxies, and to be responsible for many observational 
results including the Mbh — cr relation and the suppression of star formation in elliptical 
galaxies (e.g., Silk & Rees 1998; Fabian 1999; King 2003; Sanonov et al. 2005; Murray 
et al. 2005; Di Matteo et al. 2005; Croton et al. 2006; Ciotti & Ostriker 2001, 2007). 
While significant progresses have been made, many details of the feedback are still uncertain 
(Ostriker et al. 2010). One approach to improve this situation is to look for more direct 
observational evidence for feedback and investigate them carefully. 
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Feedback often causes intermittent activities of AGNs, thus it should be helpful to look 
for observational evidence for such an intermittency (see review by Czerny et al. 2009). 
One evidence is the central dominant galaxies in galaxy clusters observed by Chandra X- 
ray Observatory and XMM/Newton. X-ray images show cavities and ripples that gives the 
evidence for repetitive outbursts. Another example, which we will focus on in this paper, is 
radio galaxies. Many case studies of radio structure directly show evidence for two or more 
of active periods. Some sources are very compact, only a few kpc in size, which indicates 
that they are very young. These small radio sources are found to constitute 10%-30% of all 
sources in a flux-limited sample. O'Dea & Baum (1997; see also Snellen et al. 2000) have 
studied the statistical properties of a combined sample of gigahertz-peaked spectrum sources 
(GPS) and compact steep spectrum radio sources (CSS). Their results show the existence 
of far too many compact (young) sources in comparison with the number of galaxies with 
extended old radio structures. If the total activity period lasts for 10 8 yr, the number of 
sources with the ages below 10 3 yr should be roughly 10 5 times lower than the number of 
sources older than 10 3 year, which is however not the case. The most likely explanation for 
the overabundance is that the activity is intermittent, i.e., the sources undergo an active 
or outburst phase lasting for 10 4 yr which recurs every 10 5 yr (Reynolds & Begelman 1997; 
Kaiser et al. 2000; Kunert-Bajraszewska et al. 2005). Detection of several candidates for 
dying compact sources supports this view (Giroletti et al. 2005; Parma et al. 2007). 

To explain the origin of such intermittent activity, Czerny et al. (2009) have compiled 
a sample consisting of 72 GPS with measured age. The age is determined based on their 
measured hot-spot separation speed (kinematic age) or on the synchrotron cooling time. The 
majority of the age values are within the range of 200 yr and 10 4 yr. The monochromatic 
luminosity -Lsghz (= vL v at 5GHz) are all available for these sources, with typical values 
of 10 42 -10 45 erg s _1 and average value of 10 436 erg s _1 . This value can be used to estimate 
their bolometric luminosities of the accretion flows, -L a cc, i n the following way (Czerny et 
al. 2009). For some sources, the 2-10 keV luminosity L x are available. It is found that L x 
is comparable to or a factor of ~ 5 higher than, Lsghz (Vink et al. 2006; Siemiginowska 
et al. 20080 To estimate L„ c from £ we have to assume some template of the spectral 
energy distribution (SED) of the sources, either luminous AGNs or low-luminosity AGNs 
(LLAGNs), because their SEDs are characteristically different (Ho 1999). Assuming the 
SEDs of these GPS sources are similar to those of luminous AGNs, L acc is at least 10 
times higher than L x (Elvis et al. 1994). Assuming a fiducial number of 50, we have 
L acc ~ 50L X ~ 100L 5G Hz ~ 10 44 — 10 47 ergs _1 with averaged value of 3 x 10 45 ergs _1 . For 



lr The caveat is that the observed X-ray emission may not be dominated by the accretion flow (Stawarz 
et al. 2008). 
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several sources whose black hole masses are measured, L acc ~ (0.03 — 1.5)LEdd with typical 
value of ~ 0.1L Edd (Czerny et al. 2009). However, as we will argue below, these sources 
are more likely low-luminosity AGNs whose SEDs are distinctively different from luminous 
AGNs, characterized mainly by the lack of big-blue-bumps (Ho 1999; 2008). In this case, 
L acc ~ 10L X (Ho 1999), so we have L acc ~ (10 - 50)L 5 gh z ~ (10 43 - 5 x 10 46 ) ergs" 1 . The 
typical Eddington ratio is then 

La.cc ~ 0.02Z>Edd- Another statistical result from the sample 
in Czerny et al. (2009) is that the more luminous the sources are, the younger they are. 

The only model proposed so far to explain the intermittent activity of these compact 
radio sources is the thermal instability of a radiation-pressure-dominated standard thin disk 
(Czerny et al. 2009; see also Janiuk & Czerny 2011). In §2, we discuss several problems 
of this interpretation. We then propose in §3 an alternative model, i.e., a "Global Comp- 
ton scattering feedback" model, based on Yuan, Xie & Ostriker (2009). We assume that 
these radio sources are powered by hot accretion flows (advection-dominated accretion flow; 
ADAF) rather than the standard thin disk. The Compton heating to the electrons at 10 5 r s 
(r s = 2GM/c 2 is the Schwarzschild radius of the black hole with M being the black hole 
mass) by the radiation coming from the innermost region of the ADAF will cause an oscilla- 
tion of the black hole activity when L acc > 2%Lgdd- We calculate the durations of the active 
(or outburst) and inactive (or quiescent) phases and find that they are consistent with the 
observations. The last section (§4) is devoted to discussions. 



2. Criticisms on the radiation pressure instability model 

Czerny et al. (2009) assume that the accretion flow in these radio sources is described 
by the standard thin disk. We usually believe that this type of disk suffers from the thermal 
instability when the accretion rate is above a threshold, which corresponds to L ~ 0.025LEdd 
for supermassive black holes (it is ~ 0.2L E dd for stellar-mass black holes) (e.g., Svensson & 
Zdziarski 1994; Janiuk, Czerny & Siemiginowska 2002). This threshold is close to that of 
the global Compton scattering instability by coincidence. Assuming the SEDs of these radio 
sources are similar to those of luminous AGNs (refer to §1), they estimated their luminosities 
and found that all the sources in their sample are above this threshold so they concluded 
that they are thermally unstable. As a result of the instability, their accretion rates undergo 
high and low values periodically. They further assume this results in the high and low energy 
output of the jets. By adjusting parameters, this model can explain the timescales of both 
the active and inactive phases quite well. 

However, there are several questions about this model. The first one is whether jets can 
be formed from standard thin disks. To answer this question, let's look at the observations 
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of black hole X-ray binaries (BHXBs). It is believed that the physics of BHXBs and AGNs 
is the same, while the observational data in BHXBs is much better than the case of AGNs in 
terms of formation of jets from different accretion modes. An individual BHXB usually has 
five states, namely quiescent, low/hard, intermediate, high/soft, and very high (also called 
steep power-law) states. They are characterized by different spectra and timing features 
(see a review by McCintock & Remmilard 2006). Roughly speaking, the features of the soft 
state are similar to those of the luminous AGNs while the hard state is similar to LLAGNs. 
Different states are believed to be caused by different accretion modes. Specifically, the 
accretion flow in the soft state is believed to be described by the standard thin disk; while 
that in the hard state by the hot accretion flow such as an advection-dominated accretion 
flow (ADAF; Narayan & Yi 1994, 1995; for reviews on ADAFs see Narayan, Mahadevan & 
Quataert 1998 and Narayan & McClintock 2008) or luminous hot accretion flow (LHAF; 
Yuan 2001; Yuan & Zdziarski 2004). For reviews on the models of BHXBs and LLAGNs, 
readers are referred to Zdziarski & Gierlinski (2004), Done, Gierlinski & Kubota (2007), 
Narayan (2005), and Yuan (2007). Back to our question of jet formation, radio observations 
to BHXBs clearly show that jets exists in the hard state, evidenced by the flat radio spectrum 
and the elongated radio image; while both features disappear in the soft state (Fender 2006). 
This strongly indicates that jets can't be formed in the standard thin disk, but can in the 
hot accretion flows. Theoretical studies are consistent with this observational result (e.g., 
Livio, Ogilvie & Pringle 1999; Meier 2001). 

Observations find that the highest luminosity that hard states can reach is ~ (10 — 
30)%L Edd- This is roughly consistent with predictions of hot accretion flows (Yuan & 
Zdziarski 2004; Yuan et al. 2007). Above this luminosity, the cooling in the hot ac- 
cretion flow becomes so strong that the accretion flow can't remain hot and it collapses. 
On the other hand, strong radio emission has also been detected when BHXBs are more 
luminous, i.e, in the very high state. GRS 1915+105 and radio loud quasars may correspond 
to this "radio loud very high state" . But this does not imply that jets can also be formed in 
standard thin disks because the very high state can't be described by a standard thin disk. 
Actually, the accretion disk model of the very high state is still an unsolved problem. More 
importantly, detailed observations indicate that the strong radio emission comes from radio- 
emitting blobs ejected during the transition from a "hard very high state" to a "soft very 
high state", which may be modeled by a hot accretion flow-like and thin disk-like models, 
respectively (McClintock & Remillard 2006). Such kind of ejection is called "episodic jets", 
to discriminate their many different features compared to the "continuous jets" (see reviews 
in Fender & Belloni 2004; Fender, Belloni & Gallo 2004). Yuan et al. (2009) propose a mag- 
netohydrodynamical model for the formation of episodic jets, by analogy with the coronal 
mass ejection (CME) phenomenon in the Sun. Different from the "continuous jets" which 
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are formed in the presence of large-scale open magnetic fields, episodic jets are formed in the 
region of closed magnetic fields in the disk corona. According to this model, episodic jets 
can be formed in hot accretion flows and during their collapse into a cold thin disk. So both 
direct observations and the theoretical model indicate that episodic jets are not formed in 
standard thin disks. 

A question is then whether the compact radio sources correspond to the hard state or 
very high state. The model proposed in Reynolds & Begelman (1997) requires that the 
duration of the jet launching in the active phase must be ~ 10 4 yr. If the compact radio 
sources correspond to the very high state, this duration should be the timescale of the 
state transition. Unfortunately, the details of the state transition are still an open question, 
although it is believed to be associated with the collapse of the hot accretion flow to a cold 
disk. The observed timescale of the state transition ranges from a few days to ~ 100 days 
(Yu & Yan 2009; Zdziarski & Gierlinski 2004). Scaling this timescale with black hole mass 
by a factor of 10 7 gives 10 5 — 10 6 yr, which is one or two orders of magnitude longer than 
that required by Reynolds & Begelman (1997). Another point against the analogy between 
the compact radio sources and the very high state is the luminosity. As we state in §1, 
the typical luminosity of these radio sources is ~ 0.02LEdd- But on the other hand, the 
typical transition luminosity from hard to soft is significantly higher, among the range of 
~ (2 — 20)%L E dd- So it seems more likely that these compact radio sources correspond to 
the hard state rather than the very high state. But we would like to point out that it is hard 
at present to completely rule out the latter possibility. 

The second question is on the thermal instability of radiation dominated thin disk. As 
we mentioned at the beginning of this section, when the accretion rate is higher than a 
threshold so that the pressure in the accretion flow is dominated by the radiation pressure, 
the standard thin disk was thought to be thermally unstable (Shakura & Sunyaev 1976; 
Piran 1978). Time-dependent global hydrodynamical calculations (e.g., Honma et al. 1992; 
Janiuk et al. 2002; Li et al. 2007) found that the local thermal instability will result in 
the "limit-cycle" behavior. However, Gierlinsk & Done (2004) show that observations of 
the soft state of BHXBs, which has widely been believed to be described by the standard 
thin disk model, challenge the above prediction. The luminosity of the high state sources 
compiled in their sample spanes the range 0.01 < L/L Edd < 0.5, with the highest one 
well exceeding the theoretical stable limit. Thus we should expect to see the limit-cycle 
behavior. In contrast to this expectation, however, observations show little variability, which 
convincingly indicates that they are thermally stable. Recently, the non-existence of the 
thermal instability of radiation dominated thin disk has also been directly shown by three- 
dimensional MHD numerical simulation of shearing box (Hirose, Krolik & Blaes 2009). 
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Some theoretical efforts have been made to solve this puzzle. One explanation was 
proposed by Hirose, Krolik & Blaes (2009). They argued that although the stress is propor- 
tional to the total pressure^], the stress fluctuation precedes pressure fluctuations, contrary to 
the usual supposition that pressure controls the magnetic stress. This explains the thermal 
stability. An alternative explanation was proposed recently by Zheng et al. (2011). They 
argue that the previous analytical analysis of thermal stability neglect the role of magnetic 
pressure, which usually accounts for ~ 10% of the total pressure thus be dynamically im- 
portant. They show that if the magnetic pressure decreases with response to the increase of 
temperature of the accretion flow, the threshold of accretion rate above which the disk be- 
comes unstable increases significantly compared to the case of not considering the magnetic 
pressure. The physical reason is that in this case the dependence of turbulent dissipation 
heating on temperature becomes weaker. 



3. Compton scattering feedback model 

We propose that the intermittent activity of compact radio sources is caused by the 
"global Compton scattering" feedback effect in hot accretion flows. The idea of global 
Compton scattering was first proposed by Ostriker et al. (1976) and Cowie, Ostriker & 
Stark (1978) in the context of spherical accretion, and later investigated by Esin (1997), Park 
& Ostriker (2001, 2007), and Yuan, Xie & Ostriker (2009; hereafter YXO09) for rotating 
accretion flows. Here we briefly review the main results of YXO09. Most of the radiation 
comes from the innermost region of the accretion flow, ~ 10r s . The radiation will propagate 
outward and heat the electrons at large radii via Compton scattering. If the luminosity 
is high enough, the electrons will be heated above their virial temperature thus we can't 
obtain the steady solution. In this case, YXO09 argue that the black hole activity will 
oscillate between an active and inactive phases. The quantitative results are as follows. 
Given the existence of outflow in hot accretion flows, the mass accretion rate is a function 
of radius, M(r) = M out (r/r out ) s , with M out being the accretion rate at the outer boundary 
r out . In YXO09 we adopt s = 0.3, following the modeling result of the supermassive black 
hole in Sgr A* (Yuan, Quataert & Narayan 2003). We set r out = 10 5 r s , since this is the 
Bondi radius for typical parameters of black hole mass M = 1O 8 M and interstellar medium 
temperature T = 10 7 K, as in the present work. We find that when M out > M^dd, or 
equivalently L > 2%L Ed d, the electron temperature at 10 5 r s will be heated to be higher than 
the virial value thus no steady accretion solution exists. This radius is called virial radius. Its 



2 In this regard, we note that Czerny et al. (2009) assume that the stress is proportional to the geometrical 
mean between the gas and the total pressure. 
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value is found to be roughly anti-correlated with M ou t, r vir ~ 10 5 r s (M out /MEdd) _1 , mainly 
in the sense that a larger M out corresponds to a smaller r vir . Note that M out can't be too 
small, since in that case the virial radius will be out of the Bondi radius (i.e., the boundary 
of the accretion flow) thus the calculations in YXO09 may fail. 

In the present work we are satisfied with simple analytical calculations and estimations, 
but leave the detailed numerical simulation to the next step. First we note that the typical 
luminosity of the sources compiled in Czerny et al. (2009) is ~ 2%LEdd> consistent with the 
predicted threshold of luminosity in YXO09 above which the black hole activity oscillates 
due to the global Compton heating. In the following we analyze this oscillation behavior in 
more detail, calculating the timescales of the active and inactive phases and comparing with 
observations. The readers are referred to Fig. 1 for the schematic configuration of the two 
phases. 

We assume the mass of the black hole M = 10 8 MsM & . The density and temperature 
of the ISM at Bondi radius are denoted as uism and Tism ~ 10 7 K. This high temperature 
of T ISM is supported by the direct observations of some galactic centers, such as Sgr A* 
(Baganoff et al. 2003), some nearby galaxies (Pellegrini 2005), and the radio galaxy M 87 
(e.g., Di Matteo et al. 2003). Theoretically, this could be because of the shock heating due 
to the collision between stellar winds (Quataert 2004). The outer boundary of the accretion 
flow is set at the Bondi radius, r out ~ tb ~ GM/cj. ~ 10 5 r s . 

As we state above, when L ~ 2%L Edd (or equivalently M ou t ~ ^Edd), the gas at r vir ~ r B 
will be heated to be above the virial temperature thus will diffuse. The gas within r vir can 
still be accreted, which powers the high-luminosity active phase. The duration of this phase 
is determined by the accretion timescale at r vir , 

*.* ~ W'vir) ~ — ~ 3 x 10 4 M 8 yr (1) 

V r VU.l/ yMEdd J 

(Narayan, Mahadevan & Quataert 1998). For typical values of M out ~ Msdd, as indicated 
by the observed L ~ 2%LEdd, this timescale is in good agreement with the required duration 
of active phase in Reynolds & Begelman (1997). Moreover, the sample compiled in Czerny 
et al. (2009) indicates that sources with higher luminosity have younger age. This is easy 
to understand from eq. (1), because for an ADAF we roughly have L oc M„ ut thus t act ~ 
3 x 10 4 M 8 (0.1/«)(L/0.02L E dd) _1/2 yr. 

The temperature profile within r vir ~ 10 5 r s has been calculated in YXO09. Now we 
estimate the temperature of the gas beyond 10 5 r s , which is required to calculate the duration 
of the inactive phase. The temperature is mainly determined by the Compton heating, 
viscous dissipation, and bremsstrahlung cooling (we don't consider other possible heating 
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mechanisms for simplicity). Cooling is usually negligible in an ADAF. In the non-relativistic 
limit, the Compton heating reads (YXO09), 

g C omp = n e (9 x - 9 e )Fa T « n e 9 x Fcr T oc r~ 3 ' 2 . (2) 

Here 9 X = kT x /m e c 2 and T x is the radiation temperature characterizing the spectrum of 
radiation from the ADAF, F is the radiation flux, ot is Thomson cross section. The value 
of T x can be well approximated to be (YXO09), 

T x ~ 10 9 K. (3) 

Note that this value is much larger than the typical value for the spectrum of luminous 
quasar, where T x ~ 10 7 K (Sazonov et al. 2005). This is of course because of the significant 
difference between the SEDs of luminous AGNs and LLAGNs (Ho 1999; Ho 2008; Yuan, Xie 
& Ostriker 2009). The viscous heating rate is given by 

<?vis oc n e 6] /2 Hr 2 (dtt/drf oc r~ 3 ' 7 . (4) 

In both eqs. (2) and (4), the scaling of an ADAF (with outflow) n e oc r~ 3 / 2+s = r~ L2 and 
9i = kTi/rripC 2 oc r _1 are adopted. At r ~ 10 5 r s , we have g C om P ^ <2Vis, since this is why the 
electrons are heated to be above the virial temperature. Eqs. (2) and (4) then imply that 
the heating is dominated by Compton scattering at r > 10 5 r s . The evolution of electron 
temperature is then described by 

dT 

The solution is 

T e = T e0 e- At + T x (l-e~ At ). (6) 

Here T e o is the electron temperature at the beginning of the active phase and A = AlT ^m c 2 • 
The heating timescale is the duration of the active phase t act , so the spatial range that 
electrons can be heated to T x is determined by 

At act ~ 1. (7) 

This gives 



r - ( ^3 ) - 10 B r, (8) 



1/2 

\A-nm e c 2 

In the last calculation the typical values of L ~ 2%LEdd and t act from eq. (1) are adopted. 
So we obtain that the gas up to ~ 10 6 r s can be heated to Compton temperature T x after 
the active phase. Note that we neglect the bremsstrahlung cooling in our calculation, so this 
value should be regarded as an upper limit. 
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In correspondence with the increase of the temperature, the density of the gas should 
become much lower than the unheated ISM. The accretion of these low-density gas results 
in the significant decrease of mass accretion rate and the luminosity. The source then enters 
into the inactive phase. We denote the density as ni nact and estimate its value by the pressure 
balance of these gas with the surrounding ISM, whose density and temperature are denoted 
as nisM and Ti SM (~ 10 7 K): 

n ina ct = ^ISM^T" = 10 TllSM- (9) 

We estimate the mass accretion rate with the Bondi accretion theory, so 

GM . . 

r B ondi 5-, (10) 

c 

and 

M inact ~ aM Bon di ~ ^arl ondi n- mact m p c s (r Bondi ) oc T~ 3/2 n inact . (11) 

Here a is the viscous parameter, c s (rBondi) is the sound speed at rBondi- So the accretion rate 
is up to ~ 10 5 times lower in the inactive phase compared to the active phase, which is why 
the source becomes inactive. 

The Compton heated gas from r Bo ndi to 10 6 r s will cool with time. The duration of the 
inactive phase should be determined by the depletion of the heated gas due to accretion, or 
the cooling of these heated gas, whichever is shorter. The former is hard to estimate since 
10 6 r s is well out of the influence sphere of the black hole. But it must be much longer than 
the accretion timescale determined by the ADAF theory, i.e., 

tacc » ^ ~ 10 5 M 8 (-^f 1 yr. (12) 



,0.1 

The cooling timescale reads 

tcooi = = 6 x lO^VnL Yr- 

Jbrem 

Here jbrcm is the bremsstrahlung emissivity. We have T x = 10 9 K and n in act — 10~ 2 nisM) but 
what is the value of n^M? From our statement above, we know that the condition for the 
oscillation of the black hole activity is that the luminosity of the active phase is L > 2%L E dd, 
or equivalently M(10 5 r s ) > M E dd- Given that M ou t(10 5 r s ) ~ aM Bon di = 47rar Bondi c s ni S M«V 
we obtain ni SM > 10 5 (a/0.1)~ 1 M 8 " 1 cm -3 . So we have 

/ rp \ 1/2 

= 105 Mn) (lofe) yr ' (13) 

This determines the duration of the inactive phase because it is much shorter than the 
accretion timescale at 10 6 r s (eq. 12). This value is in good agreement with the duration of 
the inactive phase required in Reynolds & Begelman (1997). 
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4. Discussion: comparisons with other works 

In this paper we propose an AGN feedback mechanism occurring at r v - 1T ~ 10 5 r s from the 
central black hole, i.e., around parsec scale. The conditions for this effect to be present are 
that the accretion flow is described by hot accretion flow (such as an ADAF), and the accre- 
tion luminosity L > 2%L Edd . The influence of this effect is that the black hole will oscillate be- 
tween an active and inactive phases, which last for 3xlO 4 (O.l/a)(M/lO 8 M )(L/2%LEdd) _1 ^ 2 Y r 
and 1O 5 (q;/O.1)(M/1O 8 M0) yr ; respectively. Such kind of short-timescale oscillation should 
be common in radio sources with L > 2%L E dd- This result implies that even though the 
fueling rate of AGNs is constant at large scales, the radiation (and also matter) output from 
the AGNs oscillates. This behavior should be taken into account when we consider the effects 
of AGNs feedback at larger scales, especially those treating the small-scale accretion as the 
"sub-grid" input of their simulations. 

Now let's compare our study with other works. One is the series of works by Ciotti & 
Ostriker (2001; 2007). In their works, they focus on the radiative heating by the central AGN 
to the ISM, but on scales much larger than ours, from several pc to kpc. More physics, such 
as supernova heating, are included there. Similar to our results, they also find that when 
the luminosity of the central AGNs is high, in their case close to L^dd, oscillation occurs. 
However, the timescale of oscillation is much longer, ~ (10 7 — 10 8 ) yr. One important 
reason is that the radiation temperature adopted in their model is much smaller, T x ~ 10 7 K. 
Physically, such a low T x corresponds to the typical spectrum of a quasar (Sazanov et al. 
2005). While when the accretion rate is lower than a certain value, the accretion flow should 
make a transition from the standard thin disk to an ADAF, as argued in Narayan & Yi (1995) 
and evidenced by the transition from a soft to a hard states in BHXBs. The spectrum emitted 
by ADAFs is significantly different from that of quasars, most significantly characterized by 
the lack of big-blue-bump, as we point out in §1 (Ho 1999; Ho 2008). The corresponding 
radiation temperature of such spectra is typically T x ~ 10 9 K (YXO09). Another possibly 
important point is the inner boundary of the simulation domain. From YXO09, the global 
Compton scattering effect is already important at ~ 10 5 r s (M/M Edd ) _1 , i.e, at sub-pc or pc 
scales. But in Ciotti & Ostriker's works, they typically set the inner boundary at several 
pc. Thus some important feedback effect may be missed. Therefore, it will be interesting 
to repeat the simulations of Ciotti & Ostriker (2001; 2007), but by taking into account the 
change of T x at low M and the location of the inner boundary. This work is in progress. 

Another series of works are by Proga and his collaborators (e.g., Proga, Stone & Kallman 
2000; Proga 2007; Kurosawa & Proga 2009). Similar to our work, they focus on the sub-pc 
and pc scales. Different from our work, but similar to Ciotti & Ostriker (2001; 2007), they 
assume that the central radiation spectrum is quasar-like. Another difference from our work 
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is that they focus on cold accretion flows and higher mass accretion rates. However, since 
oscillation behavior was also found in their works, if the accretion rate decreases to such a 
low value that the accretion flow makes a transition from a thin disk to a hot accretion flow, 
T x ~ 10 9 K should again be adopted. 
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Fig. 1. — The schematic figure of the beginning of the active (a) and inactive (b) phases. At 
the active phase, the accreting gas inside 10 5 r s is virial; the gas outside this radius is heated 
by the global Compton scattering toward T x ~ 10 9 K, which is higher than the virial value. 
Thus this phase finishes when the gas inside 10 5 r s is depleted; the source then enters into 
the inactive phase. At the inactive phase, the gas is very hot, with T ~ T x ~ 10 9 K, thus 
the mass accretion rate and the luminosity is very low. This phase finishes when the hot gas 
cools down. 



